Absence of the Mott transition in infinite dimensions 
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It is shown that the half filled Hubbard model does not display the Mott metal-insulator transi- 
Q\ ' tion, in the limit of infinite dimensions. The proposed transition on the Bethe lattice is proved to 

^\ , be an artifact. In the atomic limit, the model is always in the antiferromagnetic phase, which can 

ON ' be described by the Weiss mean field model. 
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The competition between the kinetic energy and the Coulomb repulsion of electrons in a metal can cause a metal- 
insulator transition [jlj . The basic model which should describe such a transition, is the one-band Hubbard model at 
O '■ half filling, 

H = Tz £ + U T,^ - - \) (i) 

I . c=T4 
& 1 

where the hopping oc t is between (not necessarily only nearest) neighbors (ij) on a lattice of coordination number 
z, and the constant chemical potential has been incorporated in the Hamiltonian. A theoretical understanding of the 
nature of the Mott transition is still lacking. The reason is that it will occur when the gain of kinetic energy due to 
hopping oc t equals the loss of Coulomb energy oc U , where an expansion in either t or U is not permitted. However, 
it was observed that the Hubbard model simplifies but remains nontrivial in the limit of infinite coordination number 
z or dimensions oc z so that, 1/z being small already for three-dimensional lattices, it could serve as a starting 
point for an expansion in 1/z. As further evidence it was cited that the two extreme cases t — and U = Q, are 
Q ' captured exactly in the limit of infinite z || . 

In particular, it is taken for granted that the infinite dimensional Hubbard model shows a Mott transition at zero 
' temperature and a finite critical value of U, that is, a transition between a paramagnetic metal and a paramagnetic 
insulator, after frustrating the tendency to antiferromagnetic ordering. However, there is a controversy about the 
I/*-) ■ precise form of this transition Specifically, there is the question whether the Fermi liquid quasi-particle spectral 

weight vanishes at the same value of U as the Mott gap of the one-particle excitation spectrum opens. In addition 
s ! ■ the results for the quasi-particle spectral weight differ significantly already away from the proposed transition. 

It is the intention of this article, to show that the discrepancy in the results of the various approaches predicting the 
Mott transition in infinite dimensions, can be explained by the fact that they are ill-defined, and that this transition 
is absent. The proof consists of four steps: 

First it is recalled that there is no Mott gap, if the limit of infinite dimensions is taken on a Bravais lattice |jTo| , 
■ because its one-particle density of states has an infinite width. 

Therefore, secondly, a proof for the existence of the Mott transition has been suggested for the Bethe lattice at 
I i infinite z, because its one- particle density of states would be semicircular and hence would have a finite width 

In addition, it has been argued that the same semicircular density of states is realized for a lattice where the hopping 
matrix elements are distance independent, Gaussian distributed random variables P,p"2||, to fully frustrate antiferro- 
magnetic ordering. However, it will be shown that the Bethe and the fully frustrated lattice, have no well defined 
thermodynamical limit, and hence no well defined one-particle density of states. 

Thirdly, it was argued that such a semicircular density of states, or more generally a density of states of finite 
width, can be assumed ad hoc, together with a momentum independent one-particle self-energy, to simulate the limit 
of infinite dimensions P,[TT[| without referring to a particular lattice, but retaining a finite bandwidth as in finite 
dimensions. However, it will be proved here that no Hamiltonian with such properties exists, describing a limit of 
infinite z. 

Finally, it will be shown that even the qualitative features of the obtained one-particle excitation spectrum are 
incorrect, because the limit of infinite dimensions, in the paramagnetic phase, does not reproduce correctly the atomic 
limit t — ► 0. Even in the limit of infinite dimensions, the correct atomic limit of the Hubbard model ([!]) is the 
Heisenberg antiferromagnet. In particular this will invalidate the widely used second-order perturbation theory 

For the first point, consider an infinite dimensional Bravais lattice, whose one-particle density of states po never 
vanishes J14| . Because the self-energy is momentum independent for infinite dimensions, the local one-particle Green 
function G[iu> n ) = J po(e)G(iLo n ,e)de fulfills 
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G(w - *0+) = ( ^ (£) de. (2) 

J ui — e — E(w) — iO + 

If the one-particle excitation spectrum p(to) — 3G(w — iO + )/n of the local one-particle Green function G would have 
a Mott gap, G would be real valued, if cj lies in the gap. For such values of w, except for possible isolated zeroes of 
G, the self-energy E must be real as well so that 

- Po (lo - 8ffi(w)) = 0, 

which is a contradiction. Hence a Mott transition is absent on Bravais lattices in infinite dimensions. 

For the second point, in order to have a Mott transition and to preserve the finite bandwidth of finite dimensional 
lattices, usually the Bethe lattice is used, or a fully frustrated lattice where the hopping matrix elements are distance 
independent, Gaussian distributed random variables Qg. These lattices have, however, no translational invariance: 

To begin with the Bethe lattice, the reason is, that this lattice cannot be obtained as the thermodynamical limit of 
its finite counterpart, the Cayley tree of coordination number z and linear dimension N, because its surface to volume 
ratio is for large TV and z, proportional to 1 — 1/ z hence remaining finite in the limit of large N, while for a usual 
Bravais lattice this vanishes as z/N. For the Ising model, for example, this leads to a spurious phase transition of the 
Bethe-Peierls type Jl5| . 

Now, all derivations of the limit of infinite z require at some stage [[|j5] that there is a skeleton expansion of the 
self-energy Ey(o;) which becomes local oc 5ij, and hence momentum independent. For this last step translational 
invariance of the lattice is required. While for a Bravais lattice of coordination number z this holds independently of 
the boundary conditions, asymptotically in the sense, 

Ey(cj) - JyEjj(cj) = 0(1/ z) 

E«(w) - E#(w) = 0(z/N),i^j (3) 

allowing for the thermodynamical limit, the right hand side of the last equation, is of the order one for a Cayley tree. 
Hence the arguments |^,^6| leading to a momentum independent self-energy and the semicircular density of states of 
a Bethe lattice, are invalid, because it is not allowed to take the thermodynamical limit from the outset |17[] . The 
same arguments can be applied to the fully frustrated lattice, because there by construction, every lattice point is at 
the surface. 

In order to show the third point, the self-energy is supposed to be momentum independent, and a density of states 
of finite width D is assumed ad hoc to be realized by a Hamiltonian system, in infinite dimensions. The half-filled, 
paramagnetic case is considered. Then the Luttinger-Ward functional £1 describing the free energy of the system per 
lattice site, 

n = -T^2 / Po(e)[lnG _1 (iw„,e) + (iw n - e)G(iuj n ,e) 

- E ^ SM (iUn)G(iw n , e) - l]de (4) 

of the one-particle Green function G, is known (18) to have the saddle point property jqQ = if the Dyson equation 
holds. Here, denotes the z/th order self-energy, in terms of skeleton diagrams. At low temperatures, the first 
order correction to the ground state energy is therefore given by replacing G by its T = value. All but one of the 
2v Matsubara sums in the i^th order self-energy E^ can be replaced by its T = integral limits ]l9| . This can be 
performed in 2v ways, so that the last three terms in Eq. (Q) become (iui n — e — E(iu; n ))G(ia; rl , e) — 1, to cancel each 
other. The remaining term, converted into a real integral, gives the first order correction to the ground state energy 
as Cl(T) - O(0), with 

f) = — — J J dujdef(uj)po(e) arg G^;L (cl> — i0 + , e) (5) 

where arg is defined such that arg(l) = 0, and / is the Fermi function [ fl9|| . 

It the sequel, the insulating phase is assumed to exist for sufficiently large U at zero temperature, which has been 
shown ]Ti"| ] to solve Eq. (|^) for the finite width, semicircular pq. Then the one-particle excitation spectrum p(ui) has 
a gap A, around zero frequency, and the local Green function 

-*>+)=/ ^^-/e (6) 



2 



is real valued for \u>\ < A and has a Taylor expansion around zero as G(w) 
A 2 J p(e)/e 2 de of the order one. Because po has a finite width D, this and Eq. 
to and |E(w) — uj\ > D the asymptotic behavior, 



E(w) = A 2 /(aw) + 



Using the Dyson equation 



G-\u)-iO + ,e) 



E(w -i0+), 



—au>/ A 2 + . . . with positive a = 
imply for the self-energy for small 

(7) 



(8) 



then for small temperatures T <C A 2 /D, A and frequencies u> 
correction to the free energy, argG -1 is given by 



T relevant for the first order low temperature 



arg(- 



A 2 



i0+ 



) = -hQ{uj). 



Hence a negative entropy 



fi(T) - O(0) = 2Tln2 



(9) 



(10) 



results, which is a contradiction. The physical reason is of course, that there is no density of states of finite width, 
for the infinite dimensional limit of a Bravais lattice: In particular, the Mott transition on the "Bethe lattice" is an 
artifact. Any attempt to solve the infinite dimensional Hubbard model on the "Bethe lattice" will produce spurious 
results, because the "Dyson equation" does not correspond to a real Hamiltonian system, in infinite dimensions. 

As to the fourth point, it could be asked, if not at least the qualitative features of the one-particle excitation spectrum 
in the paramagnetic phase are reliable, because in the limit of infinite z, the atomic limit t — > is reproduced exactly? 
In particular it is believed [0] that two Hubbard bands around the energies ±U/2 are emerging, as U is increased. 
The numerical procedures as well as second order perturbation theory ||] suggest this. In addition, second order 
perturbation theory reproduces, at half filling, the atomic limit of the local Green function exactly B, 



G(iw n ) 1 = iuj n - (U/2) 2 /iuj r , 



E (2)( 



(11) 



so it seems that it can be used as an interpolation scheme between the two extremes t = and U = 0, and as a 
starting point for an expansion around the atomic limit. 

However, this argument contains a flaw: Eq. ( |TT| ) shows only that because the bare second order self-energy is 
exact, there is no well defined functional of the full Green functions, in terms of skeleton self-energies: The limit of 
infinite z, in the paramagnetic phase, does not reproduce correctly the atomic limit t — ► 0. The reason is that the 
ground state for t = is highly spin degenerate, so that there is no diagram technique [pOy , to justify a momentum 
independent self-energy. Even in infinite dimensions, this degeneracy must be lifted first, to arrive at the Heisenberg 
antiferromagnet |f2l| . In fact, if the Luttinger-Ward functional is expanded to second order in the skeleton self- 
energy |22]] , neither a Mott transition nor a precursor of Hubbard bands were observed for intermediate U, in the 
paramagnetic phase. This completes the proof that there is no Mott transition in infinite dimensions. 

Finally, it is focused on the infinite dimensional Heisenberg antiferromagnet as the effective Hamiltonian in the limit 
of small t/U , of the half-filled Hubbard model in infinite dimensions, with exchange coupling J cx t 2 /(U z). Its ground 
state on a hypercubic lattice is the Neel state f23| , suggesting that the infinite dimensional limit of the Heisenberg 
antiferromagnet is the Ising model, but this has not been shown explicitly so far. For a proof, the Dyson-Maleev 
representation p3] Sf = bf{l — bfb i ), — bi and S- = b^b i — 1/2 for the respective /ferromagnetic Heisenberg 
model (coupling — J) in terms of bosons is used, acting on its ground state as vacuum. The unphysical states are 
projected out by an additional, large on-site interaction cx bfbfb i b i . Now the same arguments as for the fermionic 
system p^j can be applied. In contrast to the fermionic model, the one-particle (spin wave) terms scale as 1/z, so 
that the spin-wave density of states collapses to a S function, leaving only the ferromagnetic Ising model in the limit 
of infinite z, which transforms to the antiferromagnetic Ising model if the sign of J is reversed. The Ising model itself 
reduces in infinite dimensions to the Weiss mean field theory. Corrections to order 1/z for the Heisenberg model can 
be calculated by mapping it onto an effective spin-boson model pq| . 

To conclude, what physical picture emerges from the preceding discussion? Clearly a Mott transition in the 
paramagnetic phase, as envisaged on the basis of a limit of infinite dimensions, is absent. The physical reasons are 
twofold: Firstly, although a generic finite dimensional Bravais lattice has of course a finite bandwidth, its proper 
infinite dimensional counterpart has not, thus prohibiting a Mott transition. This defect of the theory cannot be 
corrected by enforcing a finite bandwidth in infinite dimensions, because the resulting "Dyson" equation does not 
describe a meaningful Hamiltonian. 
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Secondly, there is no well denned limit of infinite dimensions for the Hubbard model in the atomic limit in the 
paramagnetic phase. The reason is, that the ground state for a vanishing kinetic energy is highly degenerate, so that 
for an expansion around the atomic limit, this intrinsic spin degeneracy has to be lifted: A paramagnetic insulating 
phase does not exist. In infinite dimensions, a simple Weiss mean field picture emerges. 

Further work on the subject will concentrate on these long range ordering processes, which in view of the arguments 
presented above, seem to dominate the Mott phenomenon in infinite dimensions. 

The author is indebted to Prof. N. Nagaosa, Dr. S. Murakami, and Dr. Chang-Ming Ho, for numerous discussions 
on the subject. The financial support by the Japan Society for the Promotion of Science is kindly acknowledged. 



N. Mott, Metal-insulator transitions (Taylor & Francis, London, 1990). 
W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989). 
A. Georges and G. Kotliar, Phys. Rev. B 45, 6479 (1992). 

A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 
R. Noack and F. Gebhardt, Phys. Rev. Lett. 82, 1915 (1999). 



R. Bulla, |sond-mat/9902290| (unpublished). 

D. Logan and P. Nozieres, Phil. Trans. R. Soc. London A 356, 249 (1998). 
S. Kehrein, Phy s. Rev. Lett. 81, 39 12 (1998). 



J. Schlipf et al, cond- mat/990226^ (unpublished). 



A. Georges and W. Krauth, Phys. Rev. Lett. 69, 1240 (1992). 
M. Rozenberg, X. Zhang, and G. Kotliar, Phys. Rev. Lett. 69, 1236 (1992). 
A. Georges and W. Krauth, Rev. Phys. B 48, 7167 (1993). 
X. Zhang, M. Rozenberg, and G. Kotliar, Phys. Rev. Lett. 70, 1666 (1993). 
E. Muller-Hartmann, Z. Phys. B 74, 507 (1989). 
T. Eggarter, Phys. Rev. B 9, 2989 (1974). 

E. Economou, Green's Functions in Quantum Physics (Springer, Berlin, 1979). 

N. Ashcroft and N. Mermin, Solid State Physics (Saunders College, Fort Worth, 1976). 

J. M. Luttinger and J. C. Ward, Phys. Rev. 118, 1417 (1960). 

J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 

A. Abrikosov, L. Gorkov, and I. Dzyaloshinski, Methods of Quantum Field Theory in Statistical Physics (Dover, New York, 
1963). 

W. Metzner, Phys. Rev. B 43, 8549 (1991). 

E. Muller-Hartmann, Z. Phys. B 76, 211 (1989). 

T. Kennedy, E. Lieb, and B. Shastry, Phys. Rev. Lett. 61, 2582 (1988). 
E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

Y. Kuramoto and N. Fukushima, J. Phys. Soc. Jpn. 67, 583 (1998). 



4 



